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Abstract 



Two factorization approaches have been proposed for single transverse spin asymmetries. One is 
the collinear factorization, another is the transverse-momentum-dependent factorization. They have 
' been previously derived in a formal way by using diagram expansion at hadron level. If the two 

factorizations hold or can be proven, they should also hold when we replace hadrons with parton 
states. We examine these two factorizations at parton level with massless partons. It is nontrivial 
to generate these asymmetries at parton level with massless partons because the asymmetries require 
f helicity-flip and nonzero absorptive parts in scattering amplitudes. By constructing suitable parton 

states with massless partons we derive the two factorizations for the asymmetry in Drell-Yan processes. 
' It is found from our results that the collinear factorization derived at parton level is not the same as 

, that derived at hadron level. Our results with massless partons confirm those derived with single 

massive parton state in our previous works. 

> 

| 1. Introduction 

Single transverse-spin asymmetries(SSA) have been observed in various experiments. An updated 
review about the phenomenology of SSA can be found in [lj. These asymmetries are expected if scat- 
tering amplitudes have nonzero absorptive parts and helicity-flip interactions are involved. The observed 
asymmetries provide a new tool to study the structure because they are sensitive to the correlations of 
partons inside a hadron and the angular momenta of these partons. In general it is difficult to predict the 
size of observed SSA in a process involved with hadrons, especially, with a transversely polarized hadron, 
because we have not enough information about the structure of hadrons. However, for SSA observed in 
processes involving large momentum transfers certain predictions can be made by using the concept of 
QCD factorization. In the factorization perturbative- and nonperturbative effects are separated. The 
nonperturbative effects are represented by various matrix elements of QCD operators and these matrix 
elements contain information about the structure of hadrons. 

Two factorization approaches have been proposed. One is the collinear factorization in which the 
partons in a hadron entering hard scattering only carry the momenta parallel or collinear to the hadron. 
Another is the transverse-momentum-dependent(TMD) factorization, where one takes transverse mo- 
menta of partons in hadrons into account. It should be noted that in the two approaches the factorization 
is derived or proposed rather formally in the sense that one works at hadron level by using the diagram 
expansion. The purpose of our work is to examine the two factorization approaches for SSA in Drell-Yan 
processes and Semi Inclusive DIS(SIDIS) at parton level with massless partons. 

The motivation of our work is based on the following: QCD factorizations, if they are proven, are 
general properties of QCD. These factorizations hold not only with hadron states but also hold when 
one replaces the hadron states with parton states. This is in the sense that the perturbatively calculable 
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parts in factorizations do not depend on hadrons and are completely determined by hard scattering of 
partons. In other words the perturbatively calculable parts or hard parts can be extracted if we replace 
hadronic states with partonic states to calculate ingredients in factorizations. This fact can be illustrated 
by the typical example of the factorization in DIS, where one can use a parton state to calculate the 
structure functions and parton distribution functions to extract the hard part. Because at the leading 
twist of the factorization only one parton is struck, a parton state with a single parton is enough for 
doing this. QCD factorizations can also be understood from the point of view of effective theories where a 
factorization is understood as a matching procedure of operators of full QCD into operators of an effective 
theory. The matching coefficients are the perturbatively calculable parts. The matching procedure is 
done with parton states, where one requires matrix elements of full QCD operators with partonic states 
to be the same of operators of an effective theory with the same partonic states multiplied with matching 
coefficients. 

There are distinct differences between two factorization approaches for SSA. In the approach of 
collinear factorization, the nonperturbative effects responsible for SSA are represented by twist-3 matrix 
elements[21 El |4j [5], or called ETQS matrix elements. These twist-3 matrix elements only contain the 
effect of spin-flip interactions. The nonzero absorptive part or T-odd effect is generated by poles of parton 
propagators in hard scattering. Applications of the collinear factorization for SSA can be found in [HE]. 
In the approach of TMD factorization, the nonperturbative effects responsible for SSA are represented by 
matrix elements containing T-odd- and spin-flip effects. These matrix elements are the so-called Sivers[8]- 
and Collins [9j functions. So far TMD factorization has been studied carefully in [lQl CEH HI CC2 [14] only 
for physical quantities which do not contain T-odd effects. In these studies the factorization is examined 
by replacing hadrons with partons and it shows that TMD parton distributions entering the factorization 
can be defined with QCD operators consistently. For SSA as T-odd effect intensive efforts has been spent 
to study how to consistently define or interpret Sivers function as a parton distribution which is gauge 
invariant and contains initial- or final state interactions[9l [T5l ll6 1 [T7 1 ll8|. Through these studies consistent 
definitions of these distributions or fragmentation functions which contain T-odd effects can be given. 
With Sivers functions SSA has been studied extensively [El EOj EH E21 [23]. The TMD factorization 
is expected to hold at low transverse momentum region. It is interesting to note that for SIDIS and 
Drell-Yan both approaches are applicable for certain kinematic region and it can be shown that the two 
approaches are equivalent [Ml (23 EB]. 

The two factorization approaches for SSA are derived rather formally in the sense that one works with 
a diagram expansion involving hadrons, where one separates diagrams involving hadrons into a parton- 
scattering part and hadronic parts. The hadronic parts are characterized by various parton density 
matrices of hadrons. With the expansion the perturbative coefficients in the two factorizations are then 
determined at the leading order of a s . As discussed in the above, the two factorizations should also 
hold when hadrons are replaced with partons, i.e., the perturbative coefficients extracted with partonic 
states should be the same as those derived with the diagram expansion at hadron level. It should be 
emphasized for a correct collinear factorization of SSA that the relevant structure function calculated 
with a parton state should be reproduced from the factorization formula by using the twist-3 matrix 
element calculated with the same parton state, even if the parton state may be simple so that one can 
not derive the complete hard part in the factorization. This is because the hard part should not depend 
on hadrons and hence it should also not depend on the parton states which is used to replace the hadrons. 
The same also holds for the TMD factorization. This will be discussed in detail when we give all results. 

We have examined the two approaches of factorization with partonic states for SSA in Drell-Yan 
processes|27] and SIDIS [28], where we have replaced the transversely polarized hadron with a quark. In 
order to have helicity-flip effect at parton level, the quark can not be taken as massless. We give the 
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quark a small mass m and calculate every element in the two factorizations to examine them. Although 
the calculation is done at leading order, but it is nontrivial because absorptive parts of amplitudes are 
needed. It is well known that SSA involved with a massive quark is proportional to the quark mass m 
[29], [30] . SSA's calculated in [27] [28] are also proportional to m. The derived perturbative coefficients 
in the factorization formulas do not depend on the quark mass m, because the twist-3 matrix element 
and the Sivers function in the corresponding factorization are also proportional to the quark mass m. 
Our results show that the formally derived TMD factorization for SSA holds at parton level, while the 
perturbative coefficient extracted from partonic results in the collinear factorization is different than that 
in the formally derived collinear factorization. In other word, the SSA calculated with the parton state can 
not be reproduced with the formally derived collinear factorization by using the twist-3 matrix element 
calculated with the same parton state. The collinear factorization can be derived with the massive parton 
state at leading but nontrivial order of a s , but the leading order result of the twist-3 matrix element 
responsible for SSA is U.V. divergent and calls for renormalization. The collinear factorization is then 
derived by interpreting the renormalization scale fx as an effective cutoff for the observed transverse 
momentum. This may be unsatisfied and may bring some problems at higher orders of a s . 

The purpose of the current work is to study the two factorizations with massless partons. We construct 
a suitable parton state of two components to replace the transversely polarized hadron. One component 
consists of one quark, another consists of one quark and one gluon. With this parton state one can 
obtain nonzero SSA and calculate nonperturbative matrix elements for extracting perturbative coefficient 
functions in the two factorizations. With massless partons we confirm all factorization results in our 
previous works [27] 128] where a single massive quark for the replacement. We establish a modified relation 
between the twist-3 matrix element and the first fci-moment of Sivers function. The relation has been 
derived before in [18, 36J, where the problem of U.V. subtractions has not been considered. 

Our work is organized as the following: In Sect. 2 we construct a partonic state which is suitable to 
study SSA at parton level with massless partons. Then we derive the leading order result of the twist-3 
matrix element with the partonic state. This result shows how the correlation between spins of partons 
generates helicity-flip effect even if the partons are massless. In Sect. 3 we give detailed results for SSA 
in Drell-Yan processes with the parton states and derive the collinear factorization, which is the same 
as derived in [27] with massive partons. In Sect. 4 we study Sivers function with the parton state and 
derive TMD factorization with our partonic results. In Sect. 5 the relation between the twist-3 matrix 
element and the first fc^-moment of Sivers function is addressed. Sect. 6 is our summary. 



2. Partonic State and Twist-3 Matrix Element 

Before we construct the suitable partonic state with massless partons to study SSA, we recall some 
standard descriptions of spin-1/2 systems from text books. We consider a system |n[A]) with the total 
spin 1/2. The system can consist of more than one partons. Its helicity is given by A = ±. Now we 
consider a forward scattering of the system through some operator O. The transition amplitude is given 
as: 

M MXl = (n[A 2 ]|0|n[Ai]>. (1) 

The system moves in the z-direction. We use Ai 2 = ± to denote the helicity of the initial- and final 
state, respectively. The transition amplitude in the helicity space is 2 x 2 matrix and can be expanded 
as: 



M 



A2A1 



a + b ■ 



A2A1 



(2) 



with <j l {i = 1, 2, 3) are Pauli matrices. The spin of the system can also be described by a spin vector 
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(s°, s) with s 2 = — 1. The two descriptions are equivalent because the transition amplitude characterized 
by the vector s can be written with the same coefficient a and the same vector b in Eq.(2): 

M(s) = \a + b- s 1 . (3) 

Single spin azimuthal asymmetries can in general appear if the non-diagonal part of M. of certain operators 

in the helicity space, i.e., ) |_ are nonzero, or, the transverse part of the vector b, i.e., b 1 and b 2 , 

nonzero. 

In this work we will use the light-cone coordinate system, in which a vector is expressed as 
= (a + ,a-,a±) = ((a + a 3 )/^, (a - a 3 )/^, a 1 , a 2 ) and a 2 ± = (a 1 ) 2 + (a 2 ) 2 . In the light-cone 
coordinate system we introduce two light-cone vectors: = (0,1,0,0) and Z M = (1,0,0,0). We define 
the totally antisymmetric tensor in the transverse space as e^" = e al3 ^ v l a np with e 0123 = 1. 
We consider the partonic state with massless partons: 

\P[X\) = |?(p,A,)[A]> +c 1 \qfa,\ q )g(k,\ g )[\]), (4) 

with the momentum specified as: 

pt+ = pv = (p+ j o, 0, 0), Vx = x P^, fc'* = (1 - x )p M = x p>*. (5) 

The state of the system is a superposition of a single quark state and a quark-gluon state. Because we 
will replace for our purpose an initial hadron with the partonic state, all partons in the state is in the 
in-state. The state has the helicity A and the total momentum P. The single quark state has the same 
helicity X q = A. For the qg-state, the total helicity is the sum A = X g + X g . The g-state and qg-state 
carries the same color index i c as given by 

\q(p,X q ))=bl(p,X q )\0), \q{p l AMK\ a ))=Tl l b)SpiA q )a\(k,X g m, (6) 

where b\ is the quark creation operator with i as the color index, is the gluon creation operator with 
a as the color index. c\ is taken as a real number. We take c\ = 1 for simplicity. Since the quark is 
massless, the non-diagonal part of the transition amplitude in Eq.(l) calculated in perturbative theory of 
QCD will obtain contributions only from the interference between the single quark- and the quark-gluon 
state: 

M+- = (q(p, +)[+}\0\q( Pl ,+)g(k, -)[-]) + (q(pi,-)g(k,+)[+]\0\q(p, -)[-]), 

M-+ = (q(p, -)[-]\0\q( Pl ,-)g(k, +)[+]) + (q(pi,+)g(k,-)[-]\0\q(p, +)[+]). (7) 

In the above the helicity of the quark is not flipped and the color of the state is summed, i.e., the color 
index i c given in Eq.(6) is summed in Eq.(7). It is clear from Eq.(7) that for SSA one needs to study the 
forward scattering of the qg-state into the q-state and the (/-state into the gg-state. 

Now we consider the twist-3 matrix element relevant to SSA. The definition with a transversely 
polarized hadron state is: 

T F ( Xl ,x 2 )e^s^ = ^ f dm^ e - W2ix2 - Xl)P+ „ mxiP+ 

2 J 47T 

• {(P,s ± \^(y 1 n) 1 + G + ^y 2 n)^(0)\P,s ± ) - (s ± - -s ± )} . (8) 
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In the above we have suppressed the gauge links along direction n between operators. These gauge 
links make the definition gauge invariant. If we replace the hadronic state with the state |-P[A]) we can 
calculate Tp perturbatively. With the definition we can find the operator for the transition in Eq.(l) as: 

0» = g s [ ^l^ e -^(^-i)^ + -^i-i^ + ^(y in ) 7 + G +M(y 2n )^( ). ( 9 ) 

J 47T 

It is straightforward to calculate the non-diagonal part of the transition amplitude at leading order of 
a s . We have: 

M%_ = ig s Tr^(N? - l)(x 2 - xi) [5(1 - x x )5(x 2 - x )e»(-) + 5(1 - x 2 )5(x x - x )e*»(+)] , 
Mt+ = ig s Tr^(N? - l)(x 2 - Xl ) [5(1 - Xl )5(x 2 - x )e»(+) + 5(1 - x 2 )S{x 1 - x )e^(-)] , (10) 

where e /i (±) is the polarization vector of the gluon. Using the representation of the polarization vector 
of the gluon moving in the z-direction: 

= --^(0, 0, l,i), = -L(0, 0, 1, -*), e*"(+) = -£*(-), (11) 

we can express the non-diagonal part with a± = (a 1 , a 2 ) in a matrix form in the 2x2 helicity space: 



M» = -gsTr^^N? - l)(x 2 - x x yio\ [5(1 - x x )5(x 2 - x ) - 5(1 - x 2 )5( Xl - x )\ , (12) 

Therefore we have the nonzero twist-3 matrix element of our partonic state at leading order of a s : 

T F (xi,x 2 ) = g s -K^J^-(N 2 - l)(x 2 - xi) [5(1 - xi)5(x 2 - x ) - 5(1 - x 2 )5(x 1 - x )\ . (13) 

This gives an example showing how the effect of helicity flip is generated by the correlation between the 
spin of a massless quark and that of a gluon. We note here at leading order Tp(x,x) = and also that 
Tp(x\, x 2 )/(x 2 — xi) is approaching to zero with x 2 approaching x\. It should be noted that there is no 
proof in QCD to show Tp(x,x) = or Tp(x,x) ^ 0. With our partonic state we have Tp(x,x) = at 
tree-level. At one-loop level Tp(x,x) becomes nonzero, as we will show. Another interesting quantity 
which will enters the collinear factorization for SSA is defined as: 

f F ( Xl ,x 2 )s1 = -i 9 4 I d ^ (r iy^-*i)r + -iy^r + 
2 J Ait 

■{(P, s ± |^(y 1 n)7 + 75 G + ^(y 2 n)^(0)|P, s ± ) - (s ± -> -s ± )} . (14) 
Replacing the hadron state with the partonic state in Eq.(4) we have the leading order result: 

7>(xi,x 2 ) = g s TT^^-(N 2 - l)(x 2 - xi) [5(1 - x 1 )5(x 2 - x ) + 5(1 - x 2 )5(x 1 - x )] . (15) 

It should be noted that our analysis in this work can be done by constructing a general partonic 
state. With such a general parton state the analysis has some similarities to these works in [3H [32] 
in different perspectives. In [31J a hadron is represented by its various light-cone components with 
hadronic wave functions. With the representation the standard Parton Distribution Function(PDF) is 
calculated in terms of these wave functions for the purpose to study the relation between scale evolutions 
of wave functions and of PDF's. If one uses the same representation to calculate the standard DIS, one 
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can extract the perturbative coefficient function of the factorization of the standard DIS. In |32j one 
represents a hadron with its various light-cone components associated with hadronic wave functions to 
study nonforward PDF's in deeply virtual Compton Scattering(DVCS). The nonforward PDF's are then 
represented in terms of hadronic wave functions. Again, one can also use the representation to calculate 
DVCS for extracting the perturbative coefficient function in the factorization of DVCS. In both cases one 
will find that the extracted perturbative coefficient functions are the same as those extracted with a single 
parton state, because in these two cases only two partons are active or entering hard scattering. Slight 
differences exist between these works of |31|, [32] and ours. In |31[ [52] the wave functions are those of a 
real hadron, while wave functions in the general parton state are arbitrary. But this makes no difference 
in extracting perturbative coefficient functions. The wave functions in |311 132] have fc^-dependence, i.e., 
the partons also carry transverse momenta. But these transverse momenta should be much smaller than 
the scale in the hard scattering and should be neglected in the hard scattering for a factorization at 
leading twist. This also implies that there will be no difference to set k± = before or after a relevant 
calculation to extract perturbative coefficient functions. 

Our study with the parton state in Eq.(4) has also certain similarity to the one-loop study of the 
structure function <?2 or gx of DIS in (33J [Ml [35] . In these studies there is an interesting problem of how 
to identify operators to define nonperturbative correlation functions. At tree-level it is straightforward 
to obtain a factorized form for e/2 in terms of a quark-quark correlation function which contains one 
bad component of the quark field. However, this factorized form can not be consistent because the bad 
components of quark fields are not dynamically independent. The bad components can be eliminated 
through the equation of motion and the factorized form can be re-expressed in terms of field correlation 
functions only with good components. In this way, 52 is actually factorized in terms of a quark- gluon 
correlation function, which is similar to Tp in Eq.(8). In our case such a problem does not exist because 
all field correlation functions for SSA do not contain any bad component of fields. After identifying the 
correct operators for 52 one can then study the one-loop correlation of 52 where one essentially has to 
study the forward Compton scattering of one parton to one parton, two to one and one to two partons. 
This is similar as we take the simple parton state in Eq.(4) to study factorizations of SSA, indicated by 
Eq.(7). 



3. SSA in Drell-Yan Processes and Its Collinear Factorization 

We consider the Drell-Yan process: 

h A {P A , s±) + h B (P B ) - 7 *(<z) + X - r +£ + + X, (16) 

where h A is a spin-1/2 polarized hadron with the transverse spin- vector s±. We take a light-cone coor- 
dinate system in which the momenta and the spin are : 

P A b = (P A> b>P A bA0), 4 = (0,0,51). (17) 

h A moves in the z-direction, i.e., Pt is the large component. The spin of h B is not observed. The 
invariant mass of the observed lepton pair is Q 2 = q 2 . The relevant hadronic tensor is defined as: 

WVV = E/ J^e^<M^,s±)M^M0)7^ (18) 



x 



and the differential cross-section is determined by the hadronic tensor as: 



da A-na 2 em Q 2 



dQ 2 d 2 q±dq+dq- 3SQ 2 



S(q 2 -Q 2 )(^-9,u)w^, S = 2PXP B . (19) 
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We are interested in the kinematical region where q\ <C Q 2 . The hadronic tensor at leading twist 
accuracy has the structure: 



-g^efs^q^W^ + (s± a e^'q'i + s± a e°i l/ q± - 9±efs± a q±pj w!p 

+q± a (eTll + £ a ±ql) q± ■ s±W® + ■■■ (20) 

In the above, we only give the tensor structures symmetric in [iv. Wj . (i = 1, 2, 3) represent T-odd effect 

related to the spin. are responsible for unpolarized cross-sections. contributes to SSA in 

the region q\ <C Q 2 which we will study. All structure functions depend on the kinematical variables 
x = q + /P^, y = q~/P B and q±. 

To study the collinear factorization of we replace the hadron hj± with the state |P[A]) in Eq. (@]) 
and the hadron Kb with an antiquark q with the momentum p 1 = (0,p~, 0, 0). The spin of q is averaged. 
Then we calculate and try to factorize it as a convolution of Tp, the antiquark distribution q(z) 

and a perturbative coefficient function TL. By the definition of one can identify the operator for the 
transition amplitude in Eq.(l). In order to have a nonzero q± the intermediate state X must contain at 
least one gluon. To generate SSA one needs to have a nonzero absorptive part. At leading order of a s a 
nonzero absorptive can only arises in the amplitude P + q — * 7* + g. The absorptive part of the amplitude 
P + q — ► 7* + g can be represented by two sets of cut diagrams, which are given in Fig.l and Fig. 3. A 
nonzero SSA is obtained through the interference of the absorptive part of the amplitude P + q — > 7* + g 
with the tree level amplitude of q + q — ► 7* + g, the later is given by the diagrams in Fig. 2. 




OP 



W ' (/) 



Figure 1: Cut diagrams for the absorptive part of the amplitude P + q — > 7* + g. The black dot is for 
the insertion of the electromagnetic current, the broken line is for the cut. 

At first look some diagrams in Fig.l are dangerous, because in these diagrams like Fig. la and Fig.lc 
one particle absorbs a collinear particle and it will give divergent contributions in the massless limit. We 
note that the diagrams in Fig.l with or without cuts must exist because of gauge invariance. It is easy to 
find that the nonzero absorptive part from cut diagrams in Fig.l only arise if the initial gluon has zero 
momentum, i.e., k + = 0. However, a massless particle state can not be defined with zero momentum. In 
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our case, a zero-momentum gluon can not be distinguished from the quark, the quark combined with the 
zero-momentum gluon should be taken as one quark state. We define our partonic state in Eq.(|4|) with 
k + 7^ 0. Therefore, we do not need to consider the absorptive parts from Fig.l. However, it is interesting 
to explore in more detail why there are no absorptive parts from Fig.l.. 



00 



(b) 



Figure 2: Diagrams for the amplitude q + q — > 7* + g. The black dot is for the insertion of the electro- 
magnetic current. 

Taking Fig. lb as an example, the absorptive part is proportional to 

= -i^-5(k + )v(p) 7 • e(A 9 ) 7 • (-P - k). (21) 

Using the polarization vector given in Eq.(ll) and v(p)j ■ p = 0, one finds that the above is proportional 
to k + 5{k + ). Hence the absorptive part of Fig. lb is zero. For Fig. la, there is some problem because one 
always has {p\ + k) 2 = 0. However, we can give to the quark a small mass and calculate the absorptive 
part first. With the massive quark, we have p± = (Px,p^,0, 0) and {p\ + k) 2 = 2k + p^ + m 2 . The 
absorptive part is then proportional to 

7 • e(Xg)u(pi) = -i—5(k+)( 1 ■ (p! + k) + m) 7 • <X g )u( Pl ). (22) 
Pi 

By using {^■p\—m)u{p\) = one finds again that the absorptive part of Fig. la is proportional to k + 5{k + ) 
and hence it is zero. One can also consider the case that the massless quark and gluon have a relative 
transverse momentum k±. The total transverse momentum of the quark and gluon is zero. In this case 
one easily finds that the cut diagram Fig. la gives no contribution. Similarly, one finds that the absorptive 
part from Fig.lc and Fig. Id are zero. However, we find that the absorptive part from Fig.le and Fig. If 
are nonzero with k + = 0. But as discussed before, we should take the partonic state with k + 7^ 0. Hence 
the absorptive part from Fig.l is always zero. In the collinear factorization derived formally, Wrp receives 
a piece of contributions as a convolution with Tp(x,x), called as soft-pole contribution, and a piece of 
contributions as a convolution of Tp(xi, X2), called hard-pole contribution. The soft pole contributions 
at parton level at the order we consider correspond actually to the contributions from Fig.l. From the 
above discussion and also the result Tp(x,x) = in the last section we can conclude that the soft pole 
contributions will not exist at leading order of a s . We will come back to this issue later. 

We now study the contributions from Fig.3. Actually, Fig. 3a, Fig. 3b and Fig.3e are also soft-pole 
contributions, because they can be nonzero only if fc~ = 0. But it is not possible to have k~ = since the 

gluon has a nonzero transverse momentum k g ± = —q± 7^ 0. So we have only three diagrams to consider. 
It is straightforward to obtain the absorptive part from these diagrams. By cutting the diagrams one 
should take care about that an extra minus sign is associated with the insertion of the electromagnetic 
current. 



v(p)j ■ e(A g )Abs 



7 • (-p - k) 
{—§ — k) 2 + ie 



Abs 



7 • (pi + k) + m 
(pi + k) 2 — m 2 + ie 
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(a) 



rrnr 



(&) 



'(d) 



(e) 




Figure 3: Cut diagrams for the absorptive part of the amplitude P + q — > 7* + 5. The black dot is for 
the insertion of the electromagnetic current, the broken line is for the cut. 



Taking the interference of Fig.3f with the complex conjugated Fig. 2a as an example, denoted as 
T^ v \fa-> which is the transition amplitude of q(p) + q(pi, X q ) + g(k, X g ) — > q(p) + q(p, X g ), we have 

Tvrfi r d^k 1 

T ^ = -w c I (#(^)^ ^ - « - we + fe - fc «) 2 > (fc - w - kg? 

■ {f abc u(p, \)TV 7 • (p - k g )Yl ■ Pl p T c j -(k g -k- p)^T a u( Pll X q ) 

■ [(k + kg)p€ a {Xg) + {-2kg + k) • C ( Xg ) Q pa + (kg - 2k) a € p{X g)}} , (23) 

where we take the spin- and color average of the initial antiquark. Here we take the opportunity to correct 
a sign error in our previous work. In |27l [28] the absorptive parts of amplitudes have been identified 
with a wrong sign. Therefore, a minus sign should be added to the results for Sivers functions and 
structure functions for SSA, where an absorptive part is involved. A minus sign should also be added to 
the corresponding formulas in the collinear factorization for SSA. In the limit q\ <C Q 2 , the contribution 
comes from when the gluon in the intermediate state is almost collinear to the quark, i.e., k^ 3> kg. 
Therefore we have for the on-shell condition for the quark crossing the cut in Fig.3f: 



5{(p + k 



1 



2p~ 



■5{k A 



A: 



k: 



2k a 



(24) 



Expanding T^ u \f a in the limit q\ <C Q 2 we obtain the leading contribution: 



fa 



I ) - — (1 - x )y/x^5(x - x )5(l - y) 

+ ... 



xp 



iq± ■ e(X g )(l + x ) - X q (l - x )e l le l (Xg)q j L 
q- = y p-. 



(25) 



The • • • in the second line represents non-leading terms in q±. From this result one can derive the 
contribution from the interference between Fig.3f and Fig. 2a to W^ v in the helicity basis as a 2 x 2 
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matrix according to Eq.(7) at leading order of a s and in the limit q\ <C Q 2 : 

Wn f a = 9T4°Ui 9 -^(N 2 - 1)^(1 - x )V2x- 6(x - xo)S(l - y) (-^) 

• [(1 + xq) — \X g \(l — xq)] + diagonal part. (26) 

In the above we have separated the contributions from Eq.(25) to the tensor into two parts. One is from 
the first term in Eq.(25), which does not depends on the spin of the quark and is proportional to a trace 
of 7-matrices without 75. Another part is from the second term and it depends on the quark spin with 
\X q \ = 1. This part is proportional to a trace of 7-matrices with 7 5 7~. We keep | X q \ to indicate the origin 

of the corresponding term. The purpose of the separation will be discussed later. One can extract 
from the above result. For other diagrams we can perform the expansion in small q± as outlined in [27]. 
We find that the contributions from other three diagrams are power-suppressed by q± in comparison with 
W^lja. Therefore we have our final result for the total 

wP(x,y,q ± ) = -^(N 2 -l)^(l-x )5(x-xo)S(l-y)[(l + x )-\X q \(l-x )] (3^) ■ (27) 

If w£ 1} can be factorized with Tp(xi,x<z) alone as assumed in [23], with our result of Tp(xi,X2) in 
the last section and the leading order result of the antiquark distribution 

q(x) = 5(1 - x) (28) 

we can derive the following factorization formula: 

W^\x,y,qj_) = 2/ S 2\2 [ — [ —q{y2)T F (x,yi)H(x/y 1 ,y/y 2 ), 

#(£i,6) = £(1-6) - 6 +0(a s ). (29) 

This verifies the result in [27], where the above factorization formula is verified with a finite quark mass 
at leading order of a s . Beside the sign problem mentioned after Eq.(23), a slight modification to the 
factorization formula in [27] is introduced here. In [27] the factor £1 = xjy\ in the numerator of the 
above H has been identified as x. This modification has no effect if we use our partonic result to test 
the factorization, because after the integration over y\ with Tp given in Eq.(13) one has y\ = 1. We 
note that the above result of Ti is different than that of the factorization formula derived formally in the 
collinear factorization approach in [23]. Using the factorization formula in [24] instead of that in Eq.(29) 
and Tp(x\,X2) in Eq.(13) with the parton state, the partonic in Eq.(27) can not be reproduced. 

It has been shown in [51 [26] that should be factorized not only with Tp(x\,X2) but also with 

Tp(xi,xz) defined in Eq.(14). These two quantities characterize different contributions corresponding to 
different projections of parton states. One involves with 7", another involves with 757". The meaning 
of projections will be explained after Eq.(31). The two different projections lead to two different contri- 
butions in One contribution is given in our in Eq.(27) with the factor (1 + x ) in [•••]. This 
contribution should be factorized with Tp. Another is given in Eq.(27) with |A g | and should be factorized 
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with Tp. Therefore we derive the following factorized form with our partonic states: 

W^(x,y,qj_) = a % „ / —[ —q{y 2 )\T F (x,yi)n{x/yi,y/y 2 )+f F (x,yi)'H(x/yi,y/y2) 

= 6(l-b) 1 + +0(a s ), H(6,e 2 )=5(l-6) + 0(a,). (30) 
I 1 - Zl)+ 

Before ending this section we briefly discuss why the soft pole contribution appears in the collinear 
factorization derived formally. In the formal approach one works directly with the hadronic tensor in 
Eq.(18) and divide diagrams contributing to the tensor into three parts: One part consists only of 
partons, and other two parts can be defined as various parton density matrices of the two initial hadrons, 
respectively. To illustrate the appearance of the soft pole contributions clearly, we take the light-cone 
gauge n ■ G = 0. In the first step of the derivation one has the multi-parton density matrix for the 
polarized hadron: 

TjifaM) = J ^^e-^-^- fcl )-^- fc M^, S ±|^(yi)^(y 2 )^(0)|P^, S± ). (31) 

The momentum k 2 — k\ is carried by the gluon entering the parton scattering. Then one makes a collinear 
expansion for the density matrix: 

■(Pa, s ± \TP(yin) 7 + Gl(y 2 n)iP(0)\P A , s ± ) + • • • , (32) 

where we have given the projection with j~ explicitly. The • • • stand for another projection with 757" 
and other irrelevant terms. It is clear from our calculation that the contribution with the projection 757" 
corresponds to the second term in Eq.(25) or that with |AJ. Assuming the above is nonzero, especially 
nonzero at kf — k 2 = 0, and combining a factor kf — k 2 from the part of parton scattering, e.g., the 
factor k + = kf — k 2 in the contribution of Fig. lb, one can write the gauge field Gj_ in the above as 
the field strength tensor G +M . At the end, one finds the soft pole contribution with kf — k 2 =0 or, 
e.g., k + = in Fig. lb. However, if we use our parton state to calculate the above density matrix, one 
finds that it is zero at kf — k 2 = 0, as indicated by Tp in Eq.(13) with x\ = x 2 - Our partonic 
at leading order of a s does not contain any soft pole contribution. This is also in the correspondence of 
Tp(x,x) = at the leading order. Therefore, with our partonic results the soft-pole contributions in the 
collinear factorization derived formally can not be identified. Beside this, Our partonic can not be 
reproduced with our partonic Tp{x\,x 2 ) by using the collinear factorization formula derived at hadron 
level. 

We have also performed calculations of the structure function for SSA in SIDIS, the same 
factorization form is obtained by replacing corresponding quantities. This confirms our previous result 
with a massive quark state in |28j . 



4. Sivers Function and TMD Factorization 

In the TMD factorization the nonperturbative quantity responsible for SSA is the Sivers function. In 
this section we perform a calculation with our partonic states to derive the TMD factorization formula 
for SSA in Drell-Yan processes. To define the Sivers function with QCD operators we introduce a gauge 
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link along the direction u with it^ = (u + , u , 0, 0): 



L u (-oo,z) 



Pexp (^—ig s J dXu ■ G(Xu + z)^J 



(33) 



The Sivers function relevant for Drell-Yan process is defined in the limit u + <C u : 

1 f dz~d 2 z± 



4 7 (2vr) 3 

■{(P,s ± \ij;(z)LU-oo,zh + L u (-oo,0)iP(0)\P,s ± ) - (s ± 



)}, (34) 



with z^ = (0, z~,z±). Beside the renormalization scale fi, the Sivers function also depends on the 
parameter Q 2 = 4(u ■ P) 2 /u 2 . With the partonic state in Eq.(|4"]) we can calculate the function. As 
discussed before, we will not consider the so called soft pole contribution because k + ^ in the partonic 
state. At leading order we have contributions from the diagrams in Fig. 4. 




(a) 



(b) 





(O 



(/) 



(c) 





/ 


/ 






i o 









(rf) 



Figure 4: The diagrams for the transition amplitude q(pi,X g ) + g(k,X g ) — > g(pi,A g ) with the operator 
for Sivers function, the broken line is for the cut. 

By taking the cut in Fig. 4a to Fig.4f one should keep in mind that the gauge link represents a particle 
moving from t = — oo and the energy flow along the gauge link is from the right to the left. Also, an 
extra minus sign associated the joining vertex of the gauge link and the quark line should be taken into 
account [27]. Although we have given possible cut for each diagram, but not every diagram can give 
nonzero contribution. E.g., the cut in Fig.4f requires in the limit u + — > that the plus component of 
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the momentum carried by the gluon in the intermediate state must be zero. This is not possible for the 
finite p±. Therefore, there is no contribution from Fig.4f at the leading order. Similarly, we find that 
the contributions from Fig. 4a, Fig. 4c, Fig.4d and Fig.4e are zero. The only nonzero contribution comes 
from the interference of Fig. 4b with Fig.4g. It is straightforward to calculate the non-diagonal part of 
the transition amplitude and to obtain Sivers function. We have: 

q±(x,p±) = ~ ^ S (yy N ^ N c ~ l)x V2x^5(x - x ). (35) 

Using the leading order result for the TMD parton distribution of the anitquark 

q(z,k ± ) = 8(l-z)5 2 (k ± ) (36) 
and the result wj 1 ' given in Eq.(27), we can derive the TMD factorization for as 



W^\zi,z 2 ,q±) = — [ d 2 k 1 _ L d 2 k 2 ± qi - 2 1J ~ q±(z 1 , k 1± )q(z 2 , k 21 _)5 2 (k l± + k 2± - q±)H, 
N c J qj_ 

H = l + 0(a s ). (37) 



This result is in agreement with that verified by using massive parton state [28J and that derived formally. 
Beyond tree-level, a soft factor should be implemented for a complete factorization. 



5. Tp{x,x) and Its Relation to Sivers Function 

Starting from the definitions of the twist-3 matrix element and Sivers function, one can formally 
derive a relation between Tp(x, x) and the first /cj_-moment of q±(x, k±) by invoking discrete symmetries. 
The relation is[T8|[36j: 

Tp(x,x) = J d 2 k_i_k 2 _q_i_(x,k±). (38) 

In deriving the relation one assumes that q_\_(x,k±) will decrease to rapidly with k± — > oo. This 
assumption implies that one discards the U.V. subtraction of the integral, i.e., the renormalization issue. 
Substituting our q±(x, k±) in Eq.(35) one sees clearly that the integral has an U.V. divergence. In our 
previous work[27| we have made an attempt to modify the relation by taking the U.V. subtraction into 
account, where we have introduced 

q±(x,b)= J cPk±e iij -- i leiq±(x,k ± ) (39) 

to make the integral convergent in the large /cj_-region. But with massive parton state we failed to establish 
a relation between q±(x,b) and Tf(x,x). From previous studies of TMD factorization for unpolarized 
cases [1011111 [i~2l 113] it can be expected in general that certain perturbative relations exist between parton 
distributions functions in collinear factorization and the corresponding parton distributions entering TMD 
factorization. The failure may indicate certain problems when one uses massive parton states to study 
problems of SSA. In this section, we will use our partonic state with massless partons to tackle the 
problem. 

We have already seen that Tp(x, x) is zero at tree-level. To obtain nonzero Tp(xi,x 2 ) at x± = x 2 = x, 
one has to consider it at the next-to-leading order. At the next-to-leading order there are too many 
diagrams because of gauge links between the operators in the definition of Tp. To reduce the number of 
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diagrams we take the light-cone gauge n ■ G = so that all gauge links becomes unit. In this gauge the 
gluon propagator reads: 

rf"n v + q u n^ s 



(40) 

+ is V n ■ q J 

In the light-cone gauge one needs in principle to give a prescription for the n • q in the denominator of 
the propagator. However, it is irrelevant here, because we will not perform any loop integration in the 
+-component. We will use the same partonic state as in Eq.(j4|) to calculate the helicity flip amplitude to 
extract Tp(x,x). In the light-cone gauge the leading contributions to Tp(x,x) are given by the diagrams 
in Fig. 5. In the representation of contributions with diagrams, we do not take the representation of cut 
diagrams. We take the representation with the product of operators in the definition of Tp as a T-order 
product. The two representations are equivalent [37] . 






(b) (c) 

t • 




(h) 



Figure 5: The leading diagrams for Tp(x,x) in the gauge n ■ G = 0. The black dots represent inserting 
field operators appearing in the definition of Tp(xi, x^)- 

As mentioned in Section 3, there will be some problem if one deals diagrams in which two collinear 
particles are combined into a divergent propagator. Here we encounter the same problem with Fig. 5c, 
Fig.5e, Fig.5f and Fig.5h. We will give the quark a small mass m to regularize the divergent propagator 
and set m = in our final result. In fact, these diagrams do not bring any trouble as we will see. 

We first consider the contribution from Fig. 5b. For Tp(x\, X2) with x\ = X2 we find that its contri- 
bution is proportional to the integral 

' ' (41) 



2tt (pi — q) 2 + ie q 2 + ie 

where q is the momentum of the gluon exchanged between the quark- and gluon line. q + is fixed as 
q + = —k + because x\ = x<i- With the fixed q + one finds that the above integral is zero. Hence 
Fig. 5b will not contribute to Tp(x,x). Because of the same reason, Fig.5d and Fig.5g will also give no 
contributions to Tp(x,x). For Fig.5e with arbitrary x\ and X2 the contribution is proportional to the 
integral 

£%t£-- 7 ' = Q. (42) 

(27r) 2 q 2 + ie {p - q) 2 + ie 

The integral is zero because of the rotation covariance in the transverse space. With the similar reason 
one finds that Fig.5c, Fig.5f and Fig.5h do not contribute to Tp(x±, X2). Therefore, only Fig.5a gives 
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nonzero contribution to T F {x,x). In the calculation we will have U.V.- and collinear divergence. We 
use dimensional regularization by taking the dimension of the transverse space as 2 — e to regularize the 
divergences. We obtain then: 



T F (x,x,n) 



N C {N 2 C - l)x V2x^5(x - x ) 



- -7 + ln4vr) +ln^ 



(43) 



In the above we have subtracted the U.V. pole and introduced the renormalization scale /i. The remaining 
pole is a collinear divergence with its corresponding scale fj, c . 
Taking our result of q±(x, k±) in the last section we obtain: 



q±(x,b) 



g s a s 



iV c (iV 2 - l)x y/2x^8(x - x ) 



7 + In 4-7T 



In 



(44) 



where the pole is for a collinear divergence. Comparing the above results, we find that T F (x,x) and 
q±(x,b) contain the same long distance effect represented by the collinear divergence. Hence a relation 
between them can be established. With our results we get: 



T F (x,x,n) 



. 2e~T 



q±(x,b) + 0(a s ). 



(45) 



This result is in certain sense of factorization, where T F (x, x) and q_\_{x, b) have the same nonperturbative 
effects represented by the collinear divergence, and the difference is a perturbative effect. It will be 
interesting to study the relation at higher orders of a s . 



6. Summary 

There are two factorization approaches for SSA. The factorization formulas in different approaches 
are formally derived with the diagram expansion at hadron level, where one separates diagrams involving 
hadrons into a parton-scattering part and hadronic parts. The hadronic parts involve hadronic states 
and are characterized by various parton density matrices of hadrons. If these factorization hold, they 
also hold when one replaces hadronic states with partonic states, because any proved factorization is 
a general property of QCD's Green functions and the perturbatively calculable parts do not depend 
on hadrons. On the other hand, it is also important to examine any factorization with partonic states 
and eventually to prove the factorization. It is a nontrivial task to examine factorizations for SSA at 
parton level, because one must have helicity-flip and nonzero absorptive parts of scattering amplitudes 
to generate SSA and the helicity of a massless quark in QCD is conserved. 

With the motivation discussed in detail in the introduction, we have made a first step to examine these 
factorization approaches for SSA in (27J EE], where we have used single massive-quark state to replace the 
transversely-polarized hadron. The nonzero mass is responsible for the needed helicity-flip. By examining 
all elements in the factorizations we have found that the formally derived TMD factorization holds at our 
parton level, while the formally derived collinear factorization has the perturbative coefficient function 
different than that derived from our partonic results. 

In this work, we construct a simple- and complicated parton states which consists massless partons 
and are suitable for studying factorization approaches for SSA. With the constructed partonic states we 
verify all factorization results in [271 [28] derived with a single massive-parton state. The results obtained 
with the simple parton state are the same as those with the complicated parton state. Because partons 
here are massless, our work gives an example how the correlation between spins of parton can change the 
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helicity of a state. With our results we have established the relation between the twist-3 matrix element 
Tp(x,x) and the first /cj_-moment of Sivers function by taking U.V. subtraction into account. 

With our results presented here, we have found that the perturbative coefficient in the collinear 
factorization of SSA is different that derived in a formal way, i.e., with the diagram expansion involving 
hadrons. Some reasons for the difference have been discussed in Sec. 3.. The difference indicates that a 
further study of the collinear factorization for SSA is needed. 

Having examined the two factorization approaches with the constructed parton states at leading but 
nontrivial order of a s , one can study higher order corrections in the two factorizations. We notice that 
recently parts of evolutions of twist-3 matrix elements have been derived [38 ] 1391 [40] . With our partonic 
state it is straightforward to derive complete evolutions of twist-3 matrix elements. It is interesting to 
note that there is a discrepancy between the result of [38, .39] and that of [40]. A part of reasons for this 
discrepancy is that the contributions to SSA and to the evolutions from the forward scattering qq — > g 
and the reversed scattering have been not considered so far in the collinear factorization. Following our 
work one can also construct a more suitable parton state to identify these contributions. A study for 
complete contributions in the collinear factorization of SSA is under the way. 
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